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We consider the production of heavy-quark pairs in the collisions of polarized and unpolarized 
on-shell photons and present, in analytic form, the fully integrated total cross sections for total 
photon spins Jz — 0, ±2 at next-to-leading-order in QCD. Phenomenological applications include 
bb production, which represents an irreducible background to standard-model intermediate-mass 
Higgs-boson production, as well as tt production. 
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I. INTRODUCTION 

It has been emphasized by many physicists that run- 
ning a future e+e" Hnear colHder (ILC) in the photon- 
photon mode is a very interesting option (see, e.g., 
Refs. [H, Q)- The high-energy on-sheU photons can be 
generated by backward Compton scattering of laser hght 
off the high-energy electron and positron bunches of the 
collider with practically no loss in energy and luminos- 
ity. In this respect, one of the most important reactions 
to consider is heavy-quark pair production in photon- 
photon collisions. A 77 collider becomes particularly im- 
portant for studies of the standard-model Higgs boson 
when its mass is below the W^W~ production threshold. 
Then, the predominant decay is — > hh. The dominant 
background to this comes from 77 — > hh, which receives 
contributions from direct and resolved photons. We leave 
aside the latter for the time being and return to this in 
Sec. |Vl The use of longitudinally polarized photons of 
equal helicity (their angular momentum being — 0) 
suppresses this background by a factor of m^/s at the 
leading order in perturbation theory Of course, 

the reason that the Jz = channel is important is that 
the Higgs signal comes entirely from it. Nevertheless, 
the above-mentioned suppression should not necessarily 
hold in general, since QCD higher-order corrections in- 
volve gluon emission, which permits the hh system to have 
Jz. 7^ 0. Therefore, the process of bottom-quark pair pro- 
duction in polarized-photon fusion would represent an ir- 
reducible background to intermediate-mass Higgs-boson 
production. Indeed, subsequent calculations of the next- 
to-leading-order (NLO) QCD corrections have confirmed 
these expectations [1, |g] . 

Furthermore, future photon colliders will become top- 
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quark factories. The data obtained there, when com- 
bined with data on top-quark production from other re- 
actions, will certainly improve our knowledge of the top- 
quark properties (see, e.g., Ref. |3]). It should also be 
noted that the NLO corrections have a large effect on 
the threshold behavior and exhibit a peculiar spin de- 
pendence in this region. 

Heavy-flavor production in photon-photon collisions 
receives contributions from direct and resolved incident 
photons. In the first case, photons behave as pointlike 
objects, interacting directly with the quarks in the hard 
scattering, while in the second case, the photon exhibits a 
complex structure involving quarks and gluons that par- 
ticipate in the hard interaction. In this paper, we present 
analytical results for the total cross sections for heavy- 
quark pair production by both polarized and unpolarized 
direct photons. The present work builds on the previous 
work of one of us 0, |3] . In Ref. [F] , differential cross sec- 
tions were calculated analytically in dimensional regular- 
ization and cast into a very compact form. We note 
that this is the only publication where complete analyt- 
ical results for polarized and unpolarized doubly differ- 
ential cross sections are presented. In Ref. jSj] , top-quark 
pair production for energies not too far above threshold 
was studied, and the fully integrated result for the so- 
called "virtual plus soft" part of the cross section was 
derived. We also note that the results presented in the 
present work constitute the Abclian part of the gluon- 
induced hadroproduction of heavy-quark pairs. 

This paper is organized as follows. Section [III explains 
our notations. In Sec. IIIH we outline our general ap- 
proach and discuss in detail our procedure and method- 
ology. In Sec. HYI we present our analytically integrated 
total cross sections. Our conclusions are summarized in 
Sec.|Vl Finally, Appendix lAlelaborates on the calculation 
of one of the most difficult double integrals. Appendix [B] 
gives expressions for the various coefficient functions that 
appear in the main text, and Appendix [Cl displays repre- 
sentations of our basis functions in terms of generalized 
Nielsen polylogarithms. 
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FIG. 1: One-loop Feynman diagrams contributing to the photon-fusion amplitude. Wavy, curly, and solid lines represent 
photons, gluons, and heavy quarks, respectively 



II. NOTATION 

For consistency, we closely follow the notations of 
Ref. 0] . The one-loop Feynman diagrams with i-channel 
topology relevant for heavy- flavor production by the scat- 
tering of two on-shell photons are depicted in Fig.[T] The 
M-channel diagrams are obtained from the depicted ones 
by crossing the incoming photon lines. Single-giuon ra- 
diation, which arises from the tree-level diagrams with a 
gluon attached in all possible ways to the heavy-quark 
line, contributes at the same order. We assign the four- 
momenta and helicities as 

lipi,Xl)+l{P2,X2)-^QiP3) + Q{Pi)+9iP5), (2.1) 

SO that pi +P2 = P3 +P4 +P5, and have pi = P2 = p^ = 
and P3 = P4 = w^, where m is the quark mass. We 
introduce the following Mandelstam variables: 

s = {pi +^'2)^ 

t = T-m" ^ {pi-p3)^-m\ (2.2) 

u = U — = {p2 — Pzf' — rr? , 

S2 = 52 - = (pi +P2- Pzf - ^ s + t + u, 

so that S2 = in the soft-gluon limit. Introducing 

v^l + -, w = ^ (2.3) 
s s -\-t 

we may write 

< = — s(l — w), u = —svw, S2 = sv{l — w). (2.4) 

In Ref. 0, the four- momentum of the gluon was in- 
tegrated out, the squared amplitudes were summed over 
the spins and colors of the final-state heavy quarks and 
averaged over the spins of the initial photons. The dif- 
ferential cross sections dAa / dvdw and da/dvdw for the 



polarized and unpolarized cases were presented analyti- 
cally, while the total cross sections Act and a were cal- 
culated numerically. In Ref. Q, these differential cross 
sections were further integrated to obtain fully analyt- 
ical result for all the terms proportional to (5(1 — w), 
i.e., those that multiply the leading-order term. How- 
ever, the hard-bremsstrahlung contributions were left 
out, for which a suitable set of parametrizations were 
constructed. It is the aim of the present work to ana- 
lytically integrate these remaining contributions, e.g. the 
expression in Eq. (30) of Ref. [51] , except for its last term, 
proportional to dai^o/dvdw. 



III. EVALUATION 

In order to obtain the analytical result for the total 
cross section, one has to perform double integrations 
over the variables v and w, as was already mentioned 
in Sec. |TT1 The explicit forms of the relevant integrals 
li {i — 1, . . . , 16) are given in Appendix C of Ref. 0]. 
However, their direct analytical evaluation turns out to 
be very complicated in general and even an unfeasible 
task in some cases. The integrals li contain logarithms 
with square roots in their arguments, and their coefficient 
functions also depend on the integration variables. In 
several cases where direct integration is possible, one ob- 
tains expressions in terms of the generalized Nielsen poly- 
logarithms [T(|]. These polylogarithms, however, contain 
long and complicated arguments that look unnatural, so 
that we decided to find some other universal representa- 
tion that would be valid for all the integrals under con- 
sideration. 

In fact, we made use of another approach to obtain the 
results. The essence of our method consists in obtaining 
the integrated result from its expansion over the variable 
it depends on, as well as in the knowledge of the ba- 
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sis functions entering the integrated result. In the past, 
such an approach was used in Ref. [ll[ for vertex- and 
propagator-type two-loop diagrams and was also applied 
to some other problems (see, e.g., Ref. [T2]). 

In our case, the result depends on the single variable 
m? / s. We find it, however, more convenient to set up 
the expansion in the heavy-quark velocity 



by using the linear ansatz 



4m^ 



(3.1) 



The procedure for obtaining the required expansions of 
the double integrals in the variable (3, by first expand- 
ing and then integrating Eq. (30) of Ref. @, was already 
discussed in detail in Ref. H and will not be addressed 
here. We only mention that, in Ref. [1], only the first 11 
terms of the expansions were obtained, which was all one 
could achieve at that time with available computer hard- 
ware resources. For our present purposes, we needed to 
greatly enlarge the depths of the expansions. Although 
this appears to be a straightforward task at first sight, it 
turned out to be a major technical challenge in practice. 
We actually needed hundreds of expansion terms to be 
able to rebuild the final integrated results. For a given 
integral, the number of expansion terms is, of course, di- 
rectly connected to the number of functions that makes 
up our basis. Thus, the main problems were, on the one 
hand, to define the smallest possible basis and, on the 
other hand, to obtain sufficiently many terms of the ex- 
pansion. Analyzing already integrated parts of the cross 
section presented in Ref. 8] and taking into account ob- 
servations made in a number of previous phenomenolog- 
ical studies, we chose our set of basis functions to be the 
complete set of harmonic polylogarithms of Remiddi and 
Vermaseren [l^. Further detailed investigation revealed, 
however, that harmonic polylogarithms alone are not suf- 
ficient, and that some nonharmonic functions should be 
added to the basis, as will be explained below. These 
functions fall into the class of multiple polylogarithms 

It is well known that Feynman amplitudes satisfy linear 
differential equations (see, e.g., Ref. IS]). In order to es- 
tablish the structure of the results, we found it very con- 
venient to consider homogeneous differential equations 
for the various integrals U , which are of the form 



fi 

E^^''^(/5)^^^(/3)=0, 



(3.2) 



where Pif' {(3) are some polynomials and k is the order of 
the homogeneous differential equation. Having typically 
150-200 coefficients of an expansion in /3, we were able 
to establish the differential equations of the above type 
for each of the li functions. As a result, we found that 
the degrees of the polynomials Pn\f3) never exceed 14 
and that the orders of the differential equations do not 
exceed 7. After having obtained the polynomials, one 
can try to solve the homogeneous differential equations 



(3.3) 



where the sum runs over all the basis functions. We re- 
mark that the first coefficients of the original /3 expan- 
sions serve as boundary conditions for our differential 
equations. Substitution of such an ansatz into Eq. (j3.2p 
leads to an algebraic system of linear equations. 

As already mentioned, not all the integrals can be given 
in terms of harmonic polylogarithms. In particular, this 
was the case for the integral Is of Ref. ^Sj , which is one 
of the most complicated ones. We explicitly integrated 
the doubly differential distribution associated with this 
function. The details are presented in Appendix A. Nev- 
ertheless, such a direct integration would be rather te- 
dious for a majority of our functions, and the changes of 
variables described in Appendix A are not universal and, 
therefore, not applicable to the other integrals. Another 
integral that cannot be expressed in terms of harmonic 
polylogarithms is Iq. 

Originally, our ansatz contained more than 100 basis 
functions. To work with such an ansatz, we needed about 
1000 expansion coefficients in the f3 series. Finally, af- 
ter some analysis, we constructed a final set of 21 basis 
functions. They are harmonic polylogarithms, except for 
three, which are discussed in Sec. IIVI With this set of 
basis functions, the number of linear equations required 
varies from several tens to a couple of hundreds, depend- 
ing on the function li considered. Typically, one needs 
about 150-200 coefficients of the /3 expansion to find the 
solution for the double integral 7^. 



IV. 



INTEGRATED RESULTS 



The unpolarized and polarized cross sections are de- 
fined in terms of crAiA2 ^-S 



(ct++ + (T+_) , 
(cr++ - cr+_) . 



(4.1) 



We parametrize the total cross section in terms of the 
polarization of the initial beams as 



] 



(4.2) 



where j — 1 — (A1A2) involves the average product of 
the photon helicities Ai and A2. According to Eq. (|4.2p . 
j — I corresponds to the unpolarized cross section dunp, 

while j — and j ^ 2 correspond to (T-(-+ and cr^ , 

respectively. 

At NLO, Eq. (142]) can be written as 



2TTa'^e*Nr r 



f(°){j,P) + CF^fW{j,P) , (4.3) 
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where eg is the fractional electric charge of the heavy 
quark Q, N^, the number of colors, and a the fine- 
structure constant. 

The Born result is well known and reads 



1-/3 
1 + 13' 
(4.4) 



The NLO result can be presented as a linear combina- 
tion of universal basis functions -F, , 



21 



(4.5) 



where all the j dependence resides in the coefficients cp-* . 

The coefficients cf ^ are given in Appendix IB] The choice 
of the basis functions Fi is not unique. We choose the 21 
basis functions Fi as follows: 

Fi = 1, F2 — Hi, F3 = -ff-i, F4 = Hi,i, 
F5 = Fq = F7 = 

Fg = Li2(/3) - Li2(-/3), Fg = HisA, 
Fio — -ff-i^i^i, Fii — 1, F12 = 

F13 — Hii^^i, Fi4 — F15 — 



Fi7 = 2 

F18 ■ 

Fi9 = 2 

Fbn = 2 



1-62 



[Li2(6)-Li2(-fe)], 



db 



In' 



1 



5 1 + 5' 



5-62 



iln2(i + 5)-iln2(l-6) 





Li2 I I - Lis 



In 2 In 



1 



6d6 



In^ 



3 + 62 2 
F21 = P [v4i(/3)-A2 (/?)]. 



In" 



(4.6) 



The H functions appearing in Eq. (|4.6p are the so-called 
harmonic polylogarithms, defined as 

Jo 1 

H±i,a,...W = J^^^^H^Ub); (4.7) 

Li2 is the dilogarithm, defined below Eq. (|4.9p : and the 
functions Ai and A2 have the following compact one-fold 
integral representations: 



Am 



^ „ ln(l-62) 1-62 

:-^[-n(l + 6) 
ln(l + (3^ + 26)] . 



We observe that only three basis functions Fig, F20, 
and F21 of Eq. (|4.6p are not expressible in terms of har- 
monic polylogarithms (|4.7|1 . We remark that the basis 
function F20 arises only from the virtual part of the cross 
section. 

We note that all the basis functions Fi of Eq. (|4.6[) can 
also be expressed via generalized Nielsen polylogarithms. 



Sn.piy) 



(n-l)!p! J t 



InP(l-ty) ln"-ii (4.9) 



with n + p = 2,3 and complicated arguments. Special 
cases include the polylogarithm of order n, Li„(?/) = 
Sn-i,i(yli and Riemann's zeta function C„ — C(n) = 
'in(l) [lol. fl6l|. We rewrite the functions Fi in terms of 



Li 



the standard generalized Nielsen polylogarithms in Ap- 
pendix[C] We note in passing that all the Fi functions can 
be expressed in terms of multiple polylogarithms 
of depth and weight 3 with simple linear arguments. 

To verify our analytical results, we compared the nu- 
merical values for the function /^^^ of Eq. (|4.3p produced 
by our Mathematica program in the polarized and un- 
polarized cases with Table 1 of Ref. [H] . There, the values 
for f^^^f are presented as functions of the variable 



2m 



(4.10) 



We found agreement on the level of better than one part 
in 10 000. Next, we compared our numbers with those 
for z = 1,2,3,4,5,10 presented in Table 1 of Ref. ^ 
dealing with the unpolarized case. The agreement was 
at the order of one part in 10 000 or better. Finally, we 
also compared our numbers with the corresponding val- 
ues for /(i)(++) and /(^'(+-) from Ref. [f. Generally, 
we were in good agreement; however, we found deviations 
for /'^'(++) by about 3% at some values of z. 

The present results form an Abelian subset of the 
non-Abelian gluon-induced NLO contributions to heavy- 
quark pair production. Recently, the total cross section 
of this subprocess was calculated analytically for unpo- 
larized gluons in Ref. [l^ using a completely different 
approach. By modifying the color structures, it is pos- 
sible to extract the unpolarized 77 cross section from 
their result. Comparing both numerically and analyti- 
cally (after expanding in /3), we find complete agreement. 
Specifically, three nonharmonic functions Fi(x), F2(x), 
and F3(x) appearing in Eqs. (13)-(15) of Ref. [l3| can 
be expressed as linear combinations of our functions Fig , 
F20, and F21. For instance, for the most complicated 
function F^lx), one has 



1-x 



15 Vl + 6a; + ; 



(4.11) 



(4.8) where X = (l-/3)/(l+/3). 



5 



V. CONCLUSIONS 

We presented, in analytic form, the integrated total 
cross sections of heavy-quark production in polarized and 
unpolarized 77 collisions at NLO in QCD. The result is 
written as a sum over bilinear products of j-dependent 
coefficient functions and j-independent basis functions, 
where j denotes the total angular momentum of the pho- 
tons. 

We checked our analytical results by reproducing, with 
great accuracy, all the numerical values listed in the rel- 
evant tables of Refs. Furthermore, we established 
agreement with the analytic NLO result for the total 
cross section of heavy-quark production via gg fusion, 
obtained just recently in Ref. by taking the Abelian 
limit. 

Using the backscattering technique, it is straightfor- 
ward to obtain polarized-photon beams of high intensity 
at the 77 option of the ILC by colliding low-energy laser 
light with polarized electron and positron beams. 

Of some concern are resolved-photon contributions. 
On the one hand, the unpolarized cross sections of 
the contributing subprocesses were already presented in 
Ref. pjj and the polarized ones may be deduced, e.g., 
from Ref. [l^. On the other hand, such contributions 
can be suppressed by operating close to the production 
threshold. In fact, we infer from Ref. [l^l that, in the 
case of &-quark production close to threshold, the resolved 
contribution only makes up a fraction of a percent of the 
full cross section. Resolved contributions may also be re- 
duced by identifying outgoing jets coUinear to one of the 
photon beams, which are a signature of resolved-photon 
events. One can also require that the energy deposited in 
the detectors be equal to the total beam energy in order 
to account for missed jets of the type mentioned above. 
From the experimental side, we are assuming only that 
heavy-quark events can be clearly identified. 

Our computer program evaluates the total cross sec- 
tions presented here in less than a second. It is publicly 
available [IJl and uses the program package HPL [2^ . 
Being implemented in Mathematica, it does not allow 
for calculations with arbitrary precision. However, with 
some additional technical modifications, arbitrary preci- 
sion could be achieved. 
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APPENDIX A: INTEGRAL 
TRANSFORMATIONS 

The contribution proportional to the integral /g of 
Ref. can be represented in the form 



h 
N 



f vdw y/x^ + VW 
dv ^=ln-^— -es{v,w), (Al) 



where iV is a known normalization constant and 
l-P 1 + 13 a 



V2 



Wi 



v{l — v) ' 



w = 1 — w, — v^w^ + 4:av{\ — vw)^ a — viV2- 

(A2) 

In the polarized case, the coefficient function eg is simply 
replaced by Aeg. The actual expressions for Aes and es 
may be found in Eqs. (B6) and (B8) of Ref. 0, respec- 
tively. 

It is convenient to change the order of integrations as 



dv / dw 



dw 



1-/32 



dv. 



(A3) 



where = 1— 4a/w. Furthermore, xs defined in Eq. (jA2 
can be represented as 



X8 = V 



ys — + 4a [ w 



so that we may substitute 



V , ^/xi + vw 
: m ■ 



1 Vm + w 



Xs — vw 



(A4) 



(A5) 
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in Eq. dU). 

Clearly, the "natural" replacement ^ — 1/v renders 
ys just linearly dependent on the new variable ^. As a 
consequence, the integral in Eq. (|Aip will be transformed 
as 



N 

R{w) = 
where eg = es/^^ 



dw R{w) 
1-/32 



^(1-t) VVs Vvs-w 



1 V ^8 + ^ - / ^ /AR^ 

In — — — —es(x,w), (A6) 
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The next step is to replace the integration variable ^ by 
the new integration variable z = so that the square 
root is removed from the logarithm. Thus, one obtains 



R{w) 



1 

2a 



, , z + w _ , , , , 

ctz In —e^yz^w)^ (A7) 



w 



where y± — 2w{l ±t — 2a) 

It is then convenient to split R{w) into the two parts 

as 



R(w) = R+iw) - R-{w), 

1 fVvT 

R±('w) — — dz \n{z ± w) es(z , w) 

2a Jjyr 



(A8) 



1 

2a 



/y+±w 



dz± In z± es{z±,w), (A9) 



where z± = z ± w, which induces a corresponding split 
of the original integral Is/N: 



h 
N 



r(+) 



(AlO) 



After exchanging the order of integrations and perform- 
ing some algebraic manipulations, we obtain 



r(±) 



— / dz± In z± / dw 68 {z± , w) 
2a Ji-H Jw+ 



-13 

(ir±ln(l±r±) / dlwefi{r±^w),{ All) 

fi Jo 



2a 



where 



z± 



w± 



1 ± r± , w± : 
2{l-r±-2a)' 



(l-r±)2 



2(1 -r± -2a)' 



(A12) 



It turns out the function es, when expressed in terms 
of the new variables r±, is greatly simplified, and so is 
the integration over the variable w. Performing the in- 
tegrals in Eq. (jAlip . most of the terms contained in eg 
yield harmonic polylogarithms H and generalized Nielsen 
polylogarithms Sa.fc- Only the most complicated terms of 
68 lead to the structures Ai{f3) and A2{l3) in Eq. (|4.8|) . 



APPENDIX B: COEFFICIENTS 

Here, we list the coefficients cf ^ appearing in Eq. (|4.5p . 
They read 

c[^^ = -19/3 + 41/33 +fi + ^(3 + 2/32 - I/33 + /3M TT^ 



- 4/3(1 + 2/3^) ln(2/32 
j 



-32/3 - ^ U + /3 + 4/32 - ^/33 ) ^2 



^-2,3 ^ 



+ 4^(7-5^1 ln2 + 24^(l.i;3 

+ ^(2 - /3 + 2/32 + 2/3^ - 2/34 -P^- 2/36)7r2 



+ 2{7 + 0^ + W- 



24 



3 + /32 



9 ± 14/3 + 7/32 - 



40 



ln2 + 8/3^(1 + /32) In /3 
, 4(3 T/?) 



5-/32 3 ^_ ^2 



1 



+ i^i^ + /3 + 8/32 - 2/33 -1- /35)^2 
8/3 



- 2 



21 + 5/32 - 



48 



48 



3 + /32 (3-h/32)2 



In 2 



4(5 + 3/32) In /3L 



'-4,7 



= Tl8 - 9/3 ± 10/32 3/3^ T 6/3'' ± 



240 



320 



5-/32 ' (5_/32)2 



48 



1 



3 + 2(3' 



± 3/3 8/32 ip 6/3^ - 8/3"^ ± 3/3^ 



8/36) In 2 + - (1 + /32 - - /3S) In /3 



+ J 
± 



^2 2^ 2 ^ 2'^ 



160 



160 
Tt 



96 



± 



96 



(5-/32)2 ^ 5-/32 ^ (3 + /32)2 3 + ^2 
-i-(32 ± 29/3 + 8^2 ± 26^3 8/3"* t 15^^) In 2 



^5,6 — 



+ -(l + /32)ln/3 
P 



±32 - 9/3 T 8/32 + 3/3^ ± U/S'' t 



240 



± 



320 



5-/32 (5 _ ^2)2 

-I- J-(8 =F 3/3 8/32 ± 6/33 - 8/34 T 3^^ - 8/3^) In 2 
2/3 



-f-(l + /32-/34-/3f^)ln/3 



160 



± 



160 



45 9 ^fl2_3 3 
^2 ^2^^ 2^ 2^^(5-/32)2^5-/32 



1 

4^ 



(32t 29/3- 



T26/33-h8/3''±15/35)ln2 



^8 

-9,16 



-10,15 



+ -(l-f /32)ln/3 

-4 - I6/32 - 12/3-* + 8j(4 + 3/32), 

21 =F i =F /3 I8/32 ± ^3 _ 7^4 _^ ^5 

P 

-t- ^(t8 - /3 t 26/32 - 34/33 _^ g^4 ^ 7^5)^ 
-18 ip i =F /3 - 24/32 ± ^3 _^ -,^Q^4 _^ ^5 
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+ -^(t8 - 5/3 t 26/3^ + 46/3^ ± 6/3^ - 13/3^), 



„0) 

-11,14 



•-12,13 



23 

± ^ ± 5/3 - 27/32 =F 5/3^ + —P^ T 5/3^ 



33 5 
3 



1 + /3 



1 / , fl-P^ 
-9^1-1, 



<Pi 



+ :5^(±20 + 12/3 ± 17/3^ + 36^^ ± /3'^ - 14/3^), 
|±|±5/3 + 21/32T5/33-y/34T5/3^ 
+ :^(±20 - 15/3 ± 17/3^ - 30/3^ ± /3'^ + 11/3^), 

z/3 



^ I +'/'i -1'^^ 



F■2^ — — 



2d 



(1) 
21 



(d~l)2^2i +^21 
1 



-17 — 



- 6/3 + 6/33 + 6/35 - i§^(l + /32), 



In (5 In p In p + Li2(5/9) — Li2 

2 \p 

-C2-2Li2(-p), 



6 
^5 



/3 



,0) 

-19 

„(i) 



-21 



-17 

3 ' 



F^i = C3-5Si,2(-p)+4Si,2 ( -- ) +2Si,2(5) 



4Si,2 +2Si 



37r In p 



20 + 30/32 _ -^Q^4 _ ^ -^4^2 _ 5^4)^ 
-^ + ^/3^-^/3^-f(15 + 6/32-5/34), 



5932(1, P) + <y52 1, 



11 



15/3+11/3^-/3^ + ^(5-/32). 



(Bl) 



Rc{si,2 Q)-Si,2(f)+2a> 

3.(-p,-9 



p(5' p 



APPENDIX C: BASIS FUNCTIONS 



5V2 ( --^,P 



1 1 

7'p 



(CI) 



As was already mentioned in Sec. IIVI all the Fi func- 
tions in Eq. (j4.6p can be written in terms of Sn.p functions 
with 71 + p = 2, 3 and some complicated arguments. The 
functions Fi with i = 1, . . . , 16 are written in terms of 
the standard harmonic polylogarithms of Remiddi and 
Vermaseren and their representations in terms of 
Nielsen polylogarithms may be found in Ref. [Tsj . The 
functions Fi (i = 17, . . . , 21) have the following forms: 



where 



\/5± 1 



3±\/5 



"7 ^2 



-lnx[Li2(/3)-Li2(-/3)]-i^i8, 
ln/31n2 a; - 3C3 + 2 [¥>i(l, x) - (pi(-l,x)] 



Fi7 
-Pi 8 

-F19 — F^ 
F« = -ln(5-62) 



The functions and (p2 are defined as 



(C2) 



(1) 



2F, 



(2) 



2F, 



(3) 
19 ' 



) = /" — lnyln(l - ay) 



2Li2(-x) + C2 + 21na;ln(l + x) 



Lis {pi) — Lis {oLx) + In a; Li2 (ax) , 



-In x 



(2) 
19 



= In 2 In 



1 + /3 
2 
1 



In 



■/5 



'/'I 



2 

1 + /3 



1-/3 



- 1-/3 
V'l I ^1 ' 



'1 ' ■ 



1-/3 



(pi(a,a; 



<P2(Q;,a;) = / — ln(l + y) ln(l + ay) 
Jx V 

= $(l,a) - $(x,aa;), (C3) 
where (see Eq. (3.15.4) of Ref. [l^) 

<^{A,B) = / ^ ln(l + v4y) ln(l + By) 
■/o V 

- Si,2(-A) + Si,2(-B) - ^ In' 4 + ^) 



F, 



(3) 
19 



In 4 



2Li2(-a;) + C2+ 3: 



- 2C3 + 8Si,2(-a;)-2(^i(-l,a;) 
+ 2(p2(-2;2", a;) - 2ip2{-z'^, x), 



F 



(2) 



+ ln^ 

B 



Si, 2 
Si, 2 



Li2 ( ^—r—] - Li2 ^ 



A 



20 - -F20 + 2^'2p 



A 
B- A 
1 + B 



Si, 2 



VA(1 + B) 

[a{i + b) 



(C4) 



In 



1-/3 
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